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ABSTRACT 



Context. In close binary systems the atmosphere of one or both components can be significantly influenced by irradiation 
from the companion. Often the irradiated atmosphere is simulated with a single-temperature approximation for the 
entire half-sphere. 

Aims. We present a scheme to take the varying irradiation angle into account by combining several separate ID models. 
This is independent of the actual code which provides the separate stellar spectra. 

Methods. We calculate the projected area of zones with given irradiation angle and use this geometrical factor to scale 
separate ID models. As an example we calculate two different irradiation scenarios with the PHOENIX code. 
Results. The scheme to calculate the projected area is applicable independent of the physical mechanism that forms these 
zones . In the case of irradiation by a primary with T=125 000 K, the secondary forms ions at different ionisation states 
for different irradiation angles. No single irradiation angle exists which provides an accurate description of the spectrum. 
We show a similar simulation for weaker irradiation, where the profile of the Ha line depends on the irradiation angle. 

Key words, radiative transfer - binaries: close - binaries: eclipsing - stars: atmospheres 



1. Introduction 

The physical conditions in a range of objects e.g. hot exoplanets and close binaries are fundamentally influenced by 
external irradiation. Often the central source is much hotter than its companion, so that the primary can be modelled as 
stand-alone, but the irradiated secondary has an atmospheric structure which differs markedly from the unirradiated case. 
At the substellar point, where the primary is closest and the irradiation is perpendicular to the secondary's surface a hot 
spot develops. For different latitudes the intensities and irradiation angles vary, leading to zones of different temperatures 
on the secondary; in an ideal case the hot spot coincides wi th the sub-stellar poin t and rings of c onstant temperatures 
form around it. Observations of cataclysmic variables (CVs) (jPavev Sz Smithlll992l) and exoplanets (jKnutson et al.|[2007l) 
indicate, that the hot spot may be displaced from the substellar point in some cases. 

An approximation for the visibil i ty of the irradiating star (as a point sour ce) on the su rface of the irradia ted one 
was al ready presented bv lEddingtonl (|1926[ ). This is described in more detail bv iKopaj ()1959l see his Sect. IV.6). iNapieil 
(jl968|) . however, has shown that these approximations are in error not only in the 'penumbral' (a partial shadow, as 
in an eclipse, between regions of complete shadow and complete illumination) regions, b ut also at full phase, which 
makes spherical treatment of the source necessary and he provides a scheme for it. Later IWoodI ([1973) presented an 
approximation for the integrat ed bolometri c flux, which is specifically designed for computer calculations, though his 
results are similar to the ones of lKopall ()1959[ ). Simulations of irradiated st ellar atmospheres started more than two decades 
ago, using o pacities for a grey at mosphere, and have improved since (iVaz fc Nordlundl fTgSS': 'Nariai & Murata' 1983; 
iNordlund fc Vaz 1990; B rett &: Sm ith 1993). These improvements i nclude metalicity effects ( Claret & Gimcncz 1992), a 



treatment of polarised light by Cranmcr (1995), non- LTE effects (jProga et al.l 119961: iKubat) 12000) and line profiles for 



non-static atmospheres (pPeraiah fc Srinivasa Raoiri998[ ) as well as limb-darkening and gravity-brightening (Alencar fc Va3 
Il997l 119991: Alencar et ahf 1999t|Claretl 20071). Nowada ys the existing models can be used to fit parameters of individual 



systems ( Wawrzvn et al.l 2009() . iBrett fc SmithI ()1993f ) introduced the idea of entropy- matching in the convections zone. 



which presumably is well-mixed and thus can transport heat horizontally. 

A treatment of irradiation is implemented in the radiation transfer code PHOENIX (jHauschildt fc Baronlll999l ) with 
the proper geometric factors, which ar e needed for extended prim aries leading to different visibilities and irradiation angles 
on any given point on the secondary () Barman et al .11200^ I2005D . As long as full 3D calculations are out of reach of the 
current computer capabilities for state of the art atmosphere modelling codes, we need a method to combine multiple ID 
simulations each calculated for an atmospheric structure appropriate for a given temperature and irradiation angle. This 
requires an integration over the inhomogeneous stellar surface of the secondary, were regions of different temperature 
are seen under a different angle. We were surprised not to find an extensive description of the geometry for close 
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binaries in the literature; only iLamzi 3 (120031^ treats a geometrically similar problem in the context of accreting young 
stars. Th e calculations are conceptionally simple, although some what ted ious. Often this task is not explicitly described 
(|Shiman skv et al. 200^ or referred to numeric methods. , Barm an et al.l (2005) present spectra for irradiated planets, 
where the dayside is modelled in ten rings. To obtain the weighting factor, they use a statistical method and distribute 
about 2000 points on the planetary surface. A surface area is than assigned to each point and the integral is performed 
as the sum over all visible points (see also ISteinacker et al.lll996[) . In contrast, we want to proceed analytically as far 
a possible. In Sect. [5] we provide ready-to- use formulae and associated IDL code to calculate the proper weighting of 
spectra calculated for different positions on the star. Using this method we simulate two exemplary spectra of irradiated 
close binaries with stellar parameters appropriate for pre-CVs with a hot sdO and analyse which new spectral features 
are revealed when the heated side is represented by a superposition of different spectra (Sect. We end with a short 
conclusion in Sect. HI 



2. The geometry 

2.1. Model assumptions 

The calculations we present in this section are purely geometrical. We calculate the projected surface area on the plane 
of the sky of a patch on a sphere, given the limiting longitude and latitude. The problem is complicated by the fact 
that we may view the sphere under an arbitrary angle and that parts may be eclipsed by the primary star, which is also 
assumed to be spherical. 

The physical situation we have in mind is a close binary system, e.g. a pre-CV or a hot Jupiter. The orbits may be 
eccentric, as long as the ephemerides or at least the separation at the time in question is known. 

The formulae in this section are independent of the physical mechanism that produces the spectral differences on the 
observed body, they have to be taken into account by the spectral synthesis code. We think of rings of constant tem- 
perature, which develop mainly due to the irradiation angle. In principle winds are expected to smooth the temperature 
gradients, but they are difficult to model in ID. In any case, the equations in this section can be used to add up different 
spectra from different patches on the star. Also, in principle, the emission can be time-depended, but in many systems 
this is not an issue. In circularised orbits with primaries of reasonably constant luminosities a stationary temperature 
distribution on the irradiated secondary will develop. If its rotation matches the orbital period, the hot spot will coincide 
with the sub-stellar point. 



2.2. The general problem 

We adopt a standard spherical coordinate system with the origin at the centre of the secondary and the positive z-axis 
intersecting the centre of the primary: 

(x\ I sin B cos (/3 \ 
y = r sin6'siniy9 (1) 
z ) \ cos9 J 

In this model lines of constant latitude 9 (measured from the z-axis) receive the same amount of incident flux from the 
primary. Points near the terminator (which lies in the z = plane) receive less incident flux than the substellar point 
because of shallower incident angles and larger distances from the primary. The problem is still rotationally symmetric 
around the z-axis, so we choose the observer to be in the a;z-plane and call the angle between the z-axis and the 
line-of-sight (which is given by the unit vector robs) w. 





/ •^obs \ 




/ smw \ 




robs = 


Vobs 


1 = 1 
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\ Zobs ) 




\ cosw / 





Figure [T] (left panel) shows a sketch of the setup. 

We normalise all distances to the radius of the irradiated secondary star (r^ = 1). The angle a between the normal 
to the surface at any point of the surface (eqn. [1] with r = 1) and the observer (eqn. [3]) is hence simply given by: 

cos a = sin cos (/5 sin a; 4- cos cos (3) 

If w is known, e.g. for binaries with measured orbital parameters, this provides a relation 93(0:) between the longitude iy9 
and the angle a for each latitude 9: 

I \ f cos a - cos 9 cos Lu\ ^ /,x 

(p(a) = arccos — — for 6* 7^ 0, w ^ (4) 

\ sm 9 sm uj J 

If the line-of-sight is parallel to the z-axis, i.e. uj ~ 0, then 9 = a for all tp. 
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Fig. 1. left: Side-view on the secondary; middle: Some latitudes are fully visible, some partially visible and some invisible. 
right: Latitudes close to grazing incidence are covered first. 



We now proceed to calculate the surface area A between the latitudes 9i and O2, which is a ring on the surface of the 
star (Fig. [1] left panel), projected on the observer's plane of the sky. We restrict the A to those parts of the ring, that 
are seen under the angles between ai and 02- 



02 i"¥'(^^Ct2) 

A{9i,02,ai,a2) — / / sin6'cosa(i^) d(/3 d0 

9i Jip(e,ai) 



(5) 



The / sin0 AO gives the surface area covered by the ring and the factor cos a accounts for the projection on the sky. 

The problem is assumed to be symmetric to the xz-plane so we integrate only over one hemisphere (a > 0) and 
multiply the result by 2. 

If a boundaries are not explicitly given, we need to take into account that only parts of the ring are on the visible 
hemisphere (Fig. [TJ middle panel). Grazing views happen under a = — ^ or +-| and we need to distinguish between full, 
partial and not visible circles of latitude. Full circles only reach values of a in the range \uj — 6\ to jcj + 0|. Therefore, the 
integration boundaries for any 'patch' to be observable are: 



\uj-9\<a< min(|w + 6l|,7r/2) 



(6) 



In the general case the integral in eqn. [S] cannot be solved analytically, but needs to be computed numerically. 

We split the integral according to the visibility of the latitude in A^ containing fully visible circles on the day-side. 



Ap containing the partially visible circles and A^ for fully visible latitudes on the night-side: 



A{ei,e2) ^ A^iOM + Ap 
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(7) 



(8) 



(9) 



(10) 



There are, however, always only full circles visible on either the day- or the nightside. The others are, together with the 
missing part of the partially visible circles, on the opposite half-sphere, turned away from the observer. 



2.3. Special case: Isotropic radiation 

If the observed spectra do not depend on the angle a, we can perform the integration over the visible if and proceed 
analytically from eqn. [5l This is the case for isotropic emission from the surface and corresponds to a situation where 
the star shows neither limb- darkening nor limb-brightening. 

For only partially visible circles the ip integration goes up to a = ^, since there is always a 'grazing shot' when the 
circle moves out of sight, simplifying eqn. 2] to: 

(^1 = arccos(— cot cot a;) (11) 
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We decompose eqn. [7] in separate integrals for full (day- or night-side) A{ and partial circles Ap. Depending on the 
values of 9i and 62 one of them or a combination of both is used according to the limits on 9 given in eqns. ISlfTOl 



U2 /"TT 

A{ = 2 ' 



(12) 



/•02 ^ 

/ / i?^ sin6'(sin0cos((9sinw + cos6'cosw)d<,9 d6' = i?^27rcosw - sin^ 
Jdi Jo 2 

Ap — 2 / i?^ sin6'(sin6'cos(^sinw + cos6'cosa;)d(/3 d0 

J 01 Jq 



= / sin^ 0sini:j2V 1 - cot^ Scot^ w + sinO cosO cosuj2['k - arccos(cot cot w)] d6 (13) 
with the upper limit on tp from eqn. 1111 The 9 integration in Eqn. p3p still has to be solved numerically. 
2.4. Coverage 

The problem is more complicated in the case of eclipsing binaries. In this case the eqns. [SlfTUl need to be extended, so 
that areas, which are covered by the primary are excluded from the integration. As one can see in Fig. [1] (right panel) 
the first latitudes to be affected are full circles seen close to grazing-incidence. 

To describe this we need a ip{9) that uses radii, distance and angle in the plane of motion to adjust the integration 
limits accordingly. 

There are two special cases, one for a primary larger than the secondary, where a full eclipse might occur during 
transit and one for a primary smaller than secondary, where the the primary only covers the center of the secondary and 
hence four integration limits are necessary to describe both visible sides during transit. 

In short, we parametrise the edge of the primary as projected on the secondary in the plane of the sky with x', y', z' 
for a circle: 



y' = - (x' + i?sina;)2 (14) 

with Tp radius of the primary and R separation between both measured from centre to centre. 

Then we transform this system with a rotation around the y-axis to our old set of coordinates x, y, z in eqn. [TJ 

(a;' \ / sin 9 cos (p cos uj + cos 9 sin oj \ 
y' I = sin sin ] (15) 

z' / \ cos 9 cos UJ — sin 9 cos tp sin uj J 

Substituting all dashed coordinates in eqn. [T3] with eqn. [TS] we get the following expression: 



sin 6* sin = sj^'p — ([cos 6' + R] sinw + sin6'cos(/3COSCL')2 (16) 

Depending whether the 9 — const ring is still a full circle, a partial circle or the primary covers the substellar point 
we get new integration limits, still taking into account eqns. ISlfTUl 

For the outer marked ring shown in the example of Fig. [1] (right panel) the new limits are (pcover,iower = |7r and 
^coverMgher = TT respectively. 

2.5. Phase and inclination 

We still need to calculate cj from the phase p and the inclination i. This can be done with the scalar product between 
the line-of-sight and the 2;-axis: 

(cospcosi\ /0\ 
sinp ^ cospsini = coscj uj = arccos(cospsini) (17) 

cos p sin i / \1 / 

2.6. IDL code 

We provide an IDL code, that implements the formulae derived in this section in elec- 
tronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or at 
http : //www . hs . uni-hamburg . de/DE/Ins/Per / Guenther / codes/geom . pro . 

3. A star in 1.5 D 

Our geometric derivation is independent of the mechanism, which causes the different zones on the star. In this section 
we provide some examples choosing a specific spectral synthesis code to illustrate the use case of the formulae above. 
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wovelength [A] wavelength [A] 



Fig. 2. left: Spectra for three different irradiation angles with irradiation by a close hot primary. The lines of C, N and 
O are labelled. (Colour in electronic version only) right: Combined spectrum for different phases and inclination i = 85° 
with irradiation by a close hot primary. (Colour in electronic version only) 



3.1. Simulation setup 

In the following we show simulated spectra of irradiated stars. The sp ectra are calculated using a modified version of 
the PHOENIX 15.04 stellar atmosphere code (jHauschildt fc Baronl[l9 99). It can co nsistently produce the incident flux 
spectrum from the hot sdOJ Aufdenbers 2001) and model the cooler MS companion (jAUard et al.ll200lh . It also includes 
an irradiation mode (.Barman et ah 2004. 2005), which treats the radiation from a primary star in the outer boundary 
condition of the radiation transport equation. H and He are treated in full NLTE. 

We combine ID models calculated for different inclination angles and compare the resulting spectra. We present two 
examples at very different temperatures. 



3.2. Limitations of the model 

PHOENIX solves the hydrodynamic equations in ID, thus, for each zone on the secondary, where we calculate a spectrum, 
the hydrodynamical structure adjusts according to the energy input from the primary star. Therefore, the temperatures, 
pressures and densities in each layer differ between the separate zones we model. This should lead to horizontal energy 
transport between the zones in the form of heat conduction or winds, but our model cannot account for that. For small 
gradients and large physical dimensions this is an acceptable approximation. It could be improved by matching the 
entropy of t h e different ring s in th e convection zone to account for horizontal mixing in this layer (jBrett'fc SmithlllQOl 
IClaretll2b0lD . iBarman et ahl (|2004D used this method for pre-CVs. However, we decided to keep the energy input at the 
lower boundary fixed according to the effective temperature of the non-irradiated secondary. Thus, the irradiated zones 
do not reach the same adiabat at depth. Entropy-matching requires a grid of models with different intrinsic effective 
temperatures for each irradiation angle, which is beyond the scope of the simple example we want to provide here. 
Because we enforce radiative equilibrium throughout the model, the total albedo is equal to 1.0 (|Eddingtonlll926[ ). 



3.3. Irradiation by an exceptionally hot primary 

First, we choose a situation which closely resembles a pre-CV composed of a hot primary sdO and a secondary MS star. 
Reviews about detached binari es, physical processe s in clo se binary systems, and g eneral three-dimensional fluid dynamics 
in binary systems are given by Claret & Gimcnez (2001) , Beer & Podsiadlowski (2002) , and 'Marsh' (200 0[). The setup i s 
similar to the UU Sge system (Abell 1966; Bond et al. 1978; PoUacco fc Bch 1993) as modelled in Wawrz vn et al.l (|2009D . 
but here we show results for an extreme situation with a much hotter primary with Tp = 125 000 K, which irradiates a 
main-sequence secondary. The centres of mass of both components are separated by 2.5 Rq, the radius of the primary 
is 0.34 Rq, and the radius of the irradiated secondary is 0.27 Rq. 

Due to the extrem irradiation conditions the models are numerically fragile. In the simulations we present the temper- 
ature correction between the iterations is at most a few %, indicating that the models are close to the physical equilibrium 
state. 

In Fig. [2] (left panel) we show spectra for three different irradiation angles (measured towards the normal of the stellar 
surface). They are all normalised to their respective local continuum. Several species can be seen in emission, because 
the atmosphere is heated from the outside. The emission lines from different ionisation stages of C, N and O are labelled. 
Conventionally, a single temperature is used as effective temperature of the heated hemisphere, but the figure shows that 
no single temperature can exist to fit a combination of these three spectra. The N II line at 4632 A is strongest for large 
irradiation angles, that is atmospheres with a comparatively low temperature, but the strong N ill lines at 4635 A and 
4642 A are most pronounced for medium irradiation angles. For oxygen the region we show includes higher states of 
ionisation. As in the case of nitrogen, the On lines are strongest in the model with an irradiation angle of 60°, whereas 
Oiv at 4633 A is strongest for the model calculated at the substellar point, where the irradiation is strongest. The same 
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Fig. 3. Ha line profiles for three different irradiation angles from a close sdO primary. (Colour in electronic version only) 



can be observed for C ill. The C ill line at 4649 A shows the same flux for the models with irradiation angle 40° and 60°, 
in contrast to the lines at 4652 A and 4653 A. The three spectra from Fig. [5] (left panel) are then combined to form a 
stellar spectrum integrated about the stellar disk according to egns. [71fT51 The zones we use cover the area 6 — 0° — 20°, 
20° — 50° and 50° — 90° (green, red and black spectrum in Fig. [2j left panel). Their effective temperatures are 41000 K, 
35 000 K and 32 000 K. As an example we combine those spectra to calculate the emission for an inclination i = 85° and 
phases and 0.3. The relative area of these three zones is 0.15, 0.46 and 0.39 for phase according to the equations 
in Sect. [2l With increasing phase the sub-stellar point rotates out of view and at phase 0.3 only the outer regions play 
a role with a relative area of 0.15 (20° — 50°) and 0.85 (50° — 90°). The contribution to the emission lines from the 
night-side is negligible because its effective temperature is only 3 400 K. Due to the irradiation the substellar point is the 
brightest part of the photosphere and when it rotates out of view, the observed luminosity of the secondary decreases. 
Figure m (right panel) shows the predicted lines normalised to the local continuum. For phase 0.3 those lines, which are 
most prominent for large values, e.g. the O ii and C ill line between 4650 A and 4654 A, are stronger than for phase 
relative to the continuum. The spectrum contains lines of all ionisation stages discussed above from N ii to C iv. Thus, 
only a model which distinguishes different zones on the day-side can predict the changes of the emission lines with the 
phase of the binary. 



3.4. Irradiation by a cooler white dwarf 

In this second example, we irradiate the secondary with the spectrum of a primary with effective temperature of only 
20 000 K. The intrinsic temperature of the secondary is 4 100 K. All other parameters are similar to the first example: 
The centres of mass of both components are separated by 2.5 Rq, the radius of the primary is 0.34 Rq, and 0.53 Rq 
for the irradiated secondary (To// = 4100 K on the night-side). We use the same 9 boundaries as above. Their effective 
temperatures are 6500 K, 6000 K and 5400 K. The lower irradiation intensity leads to better numerical stability in the 
simulations. The atmosphere of the secondary is relativley cool and does not form highly ionised ions. The most notable 
difference between the rings with different irradiation angles can be seen in the Ha line profile (Fig. [3]), which is broad 
for the sub-stellar point, where the irradiation is more intense. Again, changes in the line profile occur when the system 
rotates with respect to the observer. 

For this example we show the effect of the irradiation on the secondary's atmosphere in more detail in Fig. 21 The 
left panel shows the evolution of the temperature with the optical depth r, which is evaluated as the continuum opacity 
at 1.2^m. The atmosphere in the zone, which represents the sub-stellar point, is the hottest for each optical depth. 
The temperatures decrease with decreasing incident fiux. All three irradiated models show a temperature-inversion with 
high temperatures on the outside and a temperature minimum close to t = 0.5. In deeper layers the temperature rises 
again and then follows the unirradiated model for the same stellar parameters. In the layers before the temperature 
minimum (r = 10"'^ — 0.5) the optical depth of the models increases with increasing irradiation for a given column 
density of the atmosphere, thus the irradiation increases the opacity of the gas. This effect has been observed first by 
iBrett fc SmithI ([19931 The right panel of Fig. H describes the structure on the atmosphere in relation to the gas density. 
For comparison purposes the upper panel again shows the evolution of the temperature. For all irradiated models the 
temperature inversion occurs around Pgas = 10^ dyn cm~^. Our models are calculated to gas temperatures about one 
order of magnitude larger than that. 

The remaining two panels of Fig. 2] characterise the radiation field; The middle panel shows the internal radiation 
flux mh, int, defined as the radiation that an atmosphere with the temperature structure from the upper panel emits in 
the absence of external irradiation. The inner boundary condition matches uh, int to the blackbody radiation for the 
intrinsic, i.e., undisturbed, temperature of the secondary. 

uh, int is a net flux, thus all inward radiation, which is reflected at the inner boundary, cancels out. In contrast, the 
radiation going outwards is summed up and therefore mh. int increases monotonically outwards by about two orders of 
magnitude. The emitted energy of all layers exceeds the energy delivered from the inner boundary, so the atmosphere 
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optical depth p_J*!il 

Fig. 4. left: Temperature of optical depth for the different zones of a secondary irradiated by a cool white dwarf (see 
text for details of the model). The lines show the following situations (top to bottom): Irradiation at the substellar 
point (green), at ^ = 45° (dotted, red), ^ = 60° (dashed, black) and the night-side, i.e. no irradiation (blue). (Colour in 
electronic version only) right: The structure of the temperature and radiation field in the irradiated atmospheres. The 
upper panel shows electron temperature Tdec, the middle panel the internal radiation flux uh, int, and the lower panel 
the external radiation flux mh, ext- See text for details of the definition of the fluxes. The different zone are marked as in 
the left panel. (Colour in electronic version only) 



would cool down in the absence of an additional energy source until mh, int matches the intrinsic flux at the inner 
boundary for all layers. At low pressures uh, int is larger for irradiation at smaller fj. because these zones are hotter in 
general. The differences at the high density are due to the setup of the model. All models are calculated to a constant 
value of r, which corresponds to different gas pressures for the different models. Thus, the inner boundary, where uh. int 
is held flxed, is reached at different gas pressures. 

The lower panel shows the external radiation flux mh. oxt- Here, a full model including the internal flux and the 
irradiation is calculated. The difference between the radiation field obtained in this case and the internal radiation flux 
(uh, int) from the previous panel is uh, cxt- It is negative, because the flux is directed inwards. 

The plot does not contain reflected external irradiation, since this cancels out in the net flux ujj, oxt- The irradiation 
per unit surface area is larger for points with small fi because they are closer to the surface of the primary and the 
irradiation angle is steeper. In the outer layers of the atmosphere the optical depth is low, so only little flux is absorbed. 
Most of the heating occurs for gas pressures of 10^ — 10^ dyn cm~^. At the bottom of this zone the external radiation 
field is completely absorbed and consequently the temperature inversion occurs. For small values of /i the radiation 
penetrates slightly deeper into the atmosphere, because the path length to reach the same depth is shorter than for 
shallower incident angles. 

4. Summary and Conclusion 

In a binary system both components may infiuence each other. We set up a coordinate system, where the z-axis connects 
the centres of mass of both components and concentrate on the secondary. Assuming that the influence of the primary 
can be parameterised with the angle 9 between surface normal and z-axis we calculate the area of a ring between two 9 
values projected on the plane of the sky of an observer at arbitrary viewing position, this is meant to be used as weighting 
factor when synthesizing spectra. The formalism includes occultation of the secondary in an eclipse if appropriate. The 
pre-requisites are purely geometrical and several effects, which change the spectrum of the secondary, can in principle 
be treated according to our formalism: Irradiation, gravitation (as long as the secondary stays spherical), abundance 
anomalies, wind impact. 

We calculate synthetic spectra with the PHOENIX stellar atmosphere code for two different irradiating primaries, 
assuming negligible horizontal heat transport, so that rings of constant temperature develop around the substellar 
point. The spectra change with the irradiation angle, thus a simulation combining multiple zones is more accurate than 
representing the entire day-side with a single temperature component. We showed that simulations with multiple zones 
can be combined from established present-day stellar atmosphere models. These calculations provide a benchmark for 
future full 3D models, which should reproduce our results in the limit of negligible horizontal heat flux. 

Acknowledgements. The authors thank T. Barman for his support in setting up calculations of irradiated atmospheres with PHOENIX. H.M.G. 
acknowledges support from DLR for project number 50OR0105. ACW was supported by the DFG (Deutsche Forschungsgemeinschaft), project 
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